In order to plot this data
1) open jupyter-notebook
2) open pdos.ipynb 

Note that ithe following files made by summing files with sumpdos.x 
sumpdos.x *\(Lu\)_wfc#1\(s\) > Lu_5s_tot.dat
sumpdos.x *\(Lu\)_wfc#2\(p\) > Lu_5p_tot.dat
sumpdos.x *\(Lu\)_wfc#3\(f\) > Lu_4f_tot.dat
sumpdos.x *\(Lu\)_wfc#4\(d\) > Lu_5d_tot.dat
sumpdos.x *\(Lu\)_wfc#5\(s\) > Lu_6s_tot.dat
sumpdos.x *\(Ht\)_wfc#1\(s\) > Ht_1s_tot.dat
sumpdos.x *\(Ho\)_wfc#1\(s\) > Ho_1s_tot.dat
sumpdos.x *\(N\)_wfc#1\(s\) > N_2s_tot.dat
sumpdos.x *\(N\)_wfc#2\(p\) > N_2p_tot.dat

The tutorial on how to make this plot comes from
https://pranabdas.github.io/espresso/hands-on/pdos/
